Model independence in two dimensions and polarized cold dipolar molecules 
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We calculate the energy and wave functions of two particles confined to two spatial dimensions 
interacting via arbitrary anisotropic potentials with negative or zero net volume. The general 
rigorous analytic expressions are given in the weak coupling limit where universality or model 
independence are approached. The monopole part of anisotropic potentials is crucial in the universal 
limit. We illustrate the universality with a system of two arbitrarily polarized cold dipolar molecules 
in a bilayer. We discuss the transition to universality as function of polarization and binding energy, 
and compare analytic and numerical results obtained by the stochastic variational method. The 
universal limit is essentially reached for experimentally accessible strengths. 
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Introduction. The existence of bound quantum states 
in a potential is a problem that appears in many scientific 
disciplines. While in three dimensions it is known that 
sufficiently shallow potentials will not bind, for one- and 
two-dimensional non-positive potentials early arguments 
by Landau and Lifshitz [l[ demonstrate that at least one 
bound state is always present for arbitrarily small poten- 
tial strength. Only much later this was generalized to 
potentials of negative or zero spatial volume by Simon 
Such bound states and their structures are not only 
of formal interest but serve a great number of applica- 
tions in several areas including exotic atoms Q , excitons 
in carbon nanotubes Q or semiconducting microcavities 
Q and organic interfaces @, pairing in two-dimensional 
Fermi gases Q , localization of adatoms on surfaces [1] or 
through nanosctructuring [§] , and in population genetics 

In this letter we introduce a general formalism to cal- 
culate the leading, as well as higher, order terms for ar- 
bitrary anisotropic weak potentials. This opens the door 
for many new applications since physical systems usu- 
ally have anisotropic features that make models assuming 
cylindrical symmetry inaccurate. We test our formulas 
in the emerging field of cold polar molecules where 
layered geometries have recently been experimentally re- 
alized [12J . The long-range dipolar potential can be made 
anisotropic by tilting the external aligning field. Novel 
supcrfluid and density-wave instabilites have been pro- 
posed in this geometry [l3| . However, very little is known 
about few-body bound states in the system, particularly 
in the experimentally relevant weak coupling limit, where 
universal features emerge independent of specific models, 
and previous works have only considered the case of per- 
pendicular polarization 14H16||. We shall compare ana- 
lytic and numerical results and investigate the approach 
towards the universal limit. Wc deduce some implica- 
tions for three-body [13] as well as for many-body states 
in two dimensions (2D). 

Formulation. Consider two particles confined to two 
spatial dimensions and interacting via a pair potential, 



V(r), where r is the relative coordinate. Using polar 
coordinates, r — (x,y) — (r cos ip,r sin ip), we can write 
the Schrodingcr equation as 
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with ip the wave function, \i the reduced mass, A the 
dimensionless strength, a 2 = 2[id 2 E/h 2 , E the energy, d 
the unit of length, and s = r/d the reduced coordinate. 
The wave function can be partial-wave expanded, 
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where a m are the expansion coefficients and the 2D radial 
wave functions, f m (s), satisfy the system of equations 
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The matrix elements, V m i, 
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carry all information about the potential. For cylindrical 
potentials V m i oc 5 m i and the different m-values decouple 
in Eq. ©. 

General derivation. The regular radial solution to 
Eq. (|3|) at the origin provides the usual boundary con- 
dition for a centrifugal barrier potential, i.e. we choose 



lim s ^ s 



f m {ct,s) = 1, where wc inserted explic- 



itly the dependence on the energy parameter, a, in /,„. 
We assume that the potential, and consequently also the 
right hand side of Eq. ((3]), diverge slower than 1/s 2 when 
,s — > 0. The integral form of the equations in ((3]) is given 
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whore the boundary condition at s = is obeyed through 
the solution of the free Schrodingcr equation: 

F m (a,s) = N /iJ| TO |(as)(2/a)l ro l|m|! , (6) 

where J m is the Bessel function and where g\ m \ in Eq. ([5]) 
is the Green's function given as 



9lm\(a,s,s') = ^V^[HQ(aa)HQ(aa r ) 



(7) 



in terms of the Hankel functions, Hm ■ For bound states 
(where a = i\a\) both the free wave functions in Eq. (|6]) 
and the Green's function in Eq. ([7]) diverge at large dis- 
tances, but their combination in Eq. ([5]) must vanish ex- 
ponentially. This provides the quantization condition for 
the energy. Mathematically the condition is given as 
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This homogeneous set of linear equations © has non- 
trivial solutions only if 



det(c m ; + 5 m i) = 0, 



(10) 



where 5 m i is Kronecker's delta. 

Solution in the limit of weak binding. We solve Eq. ([5]) 
with the boundary condition ((HJ as an expansion in A 

oo oo 

frn = A "/™ ' , «m =«oE > ( U ) 

n— n— 1 

where ao is a normalization constant multiplying all a m 
for to 7^ 0. In the lowest order we find (£ from Eq. © 
using /o (0, s) = \/s I0r S M ^ 1- Then Eq. §5§ gives 
/i 0) and/W, 

/ (1) (0, S ) = -V^fos'V 00 (s')ln(s'/s)ds' , (12) 
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/ °° V ra0 («')(»') 1 - |m l-d a 
Solving Eq. (|10p in the leading order in A gives the energy 

S=-^exp(-2 7 -2/(A 2 Ao)) , (14) 

A = -J n 00 ^V o(s)f^(0,s)d S + (15) 

E^oCw^o^ds/o 00 ^^^')/^^,^)^', 

where 7 is the Eulcr's constant. Next order correc- 
tions can be determined itcratively using the calcu- 
lated am\ fm, /q 1 ^ as described in The asymptotic 



form of the radial wavefunctions for 2D and finite-range 
anisotropic potentials is given generally as 



f m (a,s) -> y/sH^ (as)S m0 , |a|s < 



1. 
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There is always a bound state for very weak potentials 
even for zero net volume, J rVdrdip < 0, and the thresh- 
old behavior of the energy is given by Eq. HH). For 
cylindrical potentials only the first term in Eq. (|15[) con- 
tributes as Vom oc Sorn, but even the leading order is still 
in general complicated. 

Dipole results. Our first application is the system of 
two polarized molecules of reduced mass (i = M/2 con- 
fined to two parallel planes separated by a distance d. 
The corresponding dipolc-dipole potential, V, projected 
to this two-dimensional geometry is 



V(r,cp) = D 2 



d 2 — 3(r cos <p cos 6* + c?sin( 
(r 2 + d 2 ) 5 /2 



(17) 



where D is the dipole moment and 9 is the polarization 
angle measured from plane. This potential is reflection 
invariant with zero net volume for any polarization. It 
has monopole, dipole, and quadrupole terms only. The 
strength is now A — > U = MD 2 /(dti 2 ), the non-zero ma- 
trix elements are V m m , V m m ±i, V m m ±2- Contributions 



to second order in U on the right hand side of Eq. (Q2 
are included by using to = 0, ±1, ±2 since higher partial 
waves only contribute to the wave function through at 
least third order in U. We here systematically restrict 
ourselves to first order, and arrive at the energy expres- 
sion 
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where the coefficients, Aq,A\ and B\, are defined by 

A = -Ml ! + isin 2 26>+ 7^ cos 4 6» , (19) 



4 c 8 32 
+0.0053 sin 2 26> cos 2 6 - 0.0033 sin 29 cos 4 9 
0.0019 cos 6 9 - M c (0.0349 sin 2 (29) (20) 
0.0054 cos 4 (6») + 0.0156A/ C cos 2 9 + 0.0343M 2 ) , 

B x = -1.204M C - cos 2 6», (21) 
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Here a is the scattering length which is a function of 
the strength and the polarization angle [i(J. The en- 
ergy very close to threshold is exponential in U 2 , as seen 
in Eq. (| 14|) . and determined by the polarization angle 
through A . The first order terms, (Ai,Bi), in U ex- 
hibit the difference in approach to threshold for the dif- 
ferent signs of the strength, U. The second order terms, 
(A2,B 2 ), are necessary to get the correct [/-independent 
pre-exponential factor in the energy. These higher or- 
der terms are sums of a large number of contributions 
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FIG. 1: The function F(U) = -U 2 \n(\E\Md 2 /(2h 2 )) for dif- 
ferent polarization angles 9. The dashed curves are calculated 
numerically and the solid lines are from Eq. flB) . 



expressed as definite integrals. The energy for 6 = ir/2 
has been considered in recent studi es II 4| — ll 61] . While our 
result is in agreement with Ref. [16[ to second order, 
the approximation in Ref. (l5l | deviates from our results 
in Eq. (|18[) in first order through the Ai-term. We give 
more details on t he g eneralization to arbitrary polariza- 
tion directions in [l9j. 

Properties of weakly-bound states. The energy ap- 
proaches zero extremely fast with vanishing strength, 
where the wave function only depends on binding energy 
through the modified Bessel function, ^(Ms). The rate 
of convergence towards these universal characteristics is 
less clear. We therefore designed a numerical method to 
investigate these structures, see [ijl] for more details. 

The reduced energy, |a 2 |, is about 10 -3 — 10 -6 when 
U ~ 1, and the numerical results in Fig. [T]are not easily 
obtained due to the exponential square dependence on 
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FIG. 2: The value (dashed) and the derivative (solid) with 
respect to U of the function, F(U) = -U 2 \n(\E\Md 2 /(2h 2 )), 
for U = as functions of polarization angle. 



U. We approach numerically the analytic straight line 
for small U, and find a very satisfying agreement for all 9. 
We note that the omitted second order correction terms 
in Eq. (|18p quickly, U ~ 0.5, become significant in order 
to reproduce the parabolic behavior found numerically. 

The binding energy is from Fig. Q] seen to increase or 
decrease with increasing 9 for negative or positive U, 
repectively. This exponential variation with 9, deter- 
mined by 2/Ao, is shown in Fig. [2] in the limit U — > 0. 
The approach to U = also varies strongly with 9 as 
illustrated in Fig. [2] by the slope which even changes sign 
at 9 0.57. 

The structure of the wave function is revealed through 
the radial components, f m , in Eq. ([2]). In the universal 
weak binding limit, U — > 0, the m = solution, K , is 
approached for all anisotropic potentials. This feature of 
universality is reached by smearing out the wave function 
over an ever increasing part of space. Outside of the 
potential, there emerges a avery small and slowly varying 
wave function whereas the contours of the potential are 
seen at smaller distances. This variation is found even for 
very small U, since the behavior near r = is essential 
to provide binding. 

The convergence of the wave function towards Kq im- 
plies that that the mean square radius in the limit is in- 
versely proportional to the energy, (r 2 ) = 2H 2 /(3M\E\), 
as shown in and illustrated in Fig. |3] The limit is 
approached from above or below depending on the poten- 
tial. A large barrier, (as for U > 0,9 = tt/2), confines the 
wave function to small distances even when the energy 
is relatively close to the threshold. However, eventually 
the wave function leaks out and the limit is approached 
from below. For U < 0,9 ps 7r/2 the barrier is absent 
and the approach for E — > is from above. Approach 
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FIG. 3: The upper part is energy times mean square radius, 
(r 2 ). The lower part is the ratio of the second moment in the 
x-direction and (r 2 ) for different polarization angles 9. The 
energy on the £-axis is obtained by varying the strength U 
around zero, shown for U < and U > as curves with and 
without points, respectively. The horizontal line at 1/3 is the 
universal limit. 



to the universal limit for E — > proves that only m = 
waves remain in analogy to halo nuclei in three dimen- 
sions. For much stronger binding E(r 2 } increases as E 
becomes large whereas (r 2 ) is confined by the radius of 
the attractive potential. 

Another feature of universality is that the deviation 
from cylindrical symmetry on average disappears when 
E — > as shown in Fig. [31 The rate of disappearance 
varies strongly with polarization angle. For 9 = n/2 the 
potential itself is already cylindrical, whereas the largest 
deformation of potential and wave function occurs for 9 = 
0. Increasing the angle results in less deformation for all 
energies. The elongation is largest along the polarization 
direction for all U > 0. The deformations for U < are 
much smaller for comparable energies. 

Perspectives. The two-body system has to be thor- 
oughly understood to provide a solid ground for extrac- 
tion of properties of the iV-body system. We present an 
alternative proof that an arbitrary anisotropic 2D poten- 
tial with zero net volume has at least one bound state 
for infinitesimal potential strength and rigorously derive 
a general expression for the energy for very weak poten- 
tials where the system reaches universality with a wave 
function given as a modified Besscl function entirely de- 



termined by the binding energy. The analytic result for 
the energy is worked out in details for dipolar molecules, 
and we illustrate the result by accurate numerical compu- 
tations. High-order terms in the potential strength are 
necessary to get accurate values. We investigate char- 
acteristics of universality, that is approach to cylindrical 
symmetry, where the monopole component dominate and 
the mean square radius becomes inversely proportional to 
the binding energy. The universal limit is reached when 
the reduced strength, \U\, decreases below unity, which 
is the regime of current experiments [l2| . 

These results imply that three particles in 2D without 
quantum statistical restrictions would have at least two 
bound states for weak couplings with universal features 
[20l | . In conclusion, reliable predictions of many-body 
structures, including regions of various phases, must take 
these few-body structures into account as they persist to 
arbitrarily small potential strength. 
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